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Experiment  and  calculations  are  used  to  show  that  the  gas-phase  acidity  of  uracil  is  comparable  to  that  of 
HCl.  The  gas-phase  acidity  of  uracil  (denoted  here  by  U)  was  bracketed  by  proton-transfer  measurements 
involving  U  and  various  reference  acids  (denoted  here  by  A)  of  known  gas-phase  acidity.  Rate  constants  for 
proton  transfer  from  the  reference  acid  A  to  the  conjugate  anion  of  uracil,  (U-H)”,  were  measured  in  a  selected- 
ion  flow  tube  at  298  K.  Rate  constants  for  proton  transfer  from  U  to  ions  (A-H)“  were  measured  at  467  K 
in  a  flowing-afterglow  Langmuir  probe  apparatus.  Here,  (U-H)  or  (A-H)  indicates  a  U  or  A  molecule  which 
is  missing  an  H  atom,  respectively.  The  result  is  A//°acid(uracil)  =  333  =b  5  kcal  mol’^  and  AG°acid(uracil)  = 
326  i  5  kcal  mol  ^  at  298  K,  which  agrees  with  earlier  work.  Thermal  electron  attachment  to  uracil  was 
found  to  be  too  slow  to  permit  measurement  of  a  rate  constant,  consistent  with  the  gas-phase  acidity  given 
above.  G3  and  G3(MP2)  calculations  are  reported  for  uracil,  and  for  the  each  of  the  (U-H)  radicals  and 
(U-H)"  ions  that  result  from  H  or  H"^  loss  from  each  of  the  four  hydrogen  sites  of  U  (on  the  Nl,  N3,  C5,  and 
C6  positions).  From  the  calculated  total  energies  we  obtain  the  gas-phase  acidity  of  uracil,  the  four  U-H 
homolytic  bond  strengths,  and  the  electron  affinities  of  the  four  possible  fragment  radicals.  We  confirm  earlier 
work  that  the  most  acidic  site  in  uracil  is  at  the  Nl  site;  this  site  is  where  uracil  becomes  covalently  bonded 
to  a  carbon  of  the  ribose  sugar  in  RNA.  G3  calculations  for  the  Nl  site  at  298  K  give  A//°acid(uracil)  =  334.5 
kcal  mol  ^  and  AG°acid(uracil)  =  327.1  kcal  mol”^  at  298  K,  in  good  agreement  with  the  experiment.  The 
weakest  H-atom  bond  enthalpy  (at  the  Nl  site)  is  calculated  to  be  101.8  kcal  mol"^ 


I.  Introduction 

Uracil  is  one  of  the  five  nucleobases  that  pair  up  in  nucleic 
acids.  Uracil  has  acidic  character  as  well,  and  indeed  the 
hydrogen  bonding  between  uracil  and  adenine  in  nucleic  acids 
is  a  consequence  of  the  acid/base  character  of  these  important 
biological  compounds.  We  will  denote  the  uracil  molecule  by 
U.  Proton  loss  from  uracil  leaves  behind  a  negative  ion  which 
we  will  denote  by  (U-H)".  Several  experimental and 
theoretical  studies^ have  been  published  in  recent  years  on 
negative  ions  formed  from  gas-phase  uracil  and  other  nucleotide 
bases,  and  in  a  few  cases  the  acidities  were  a  focus.^"^^  Most 
notable  is  the  elegant  experiment  of  Kurinovich  and  Lee^  in 
which  the  gas-phase  acidity  of  U  was  measured  for  proton  loss 
from  the  Nl  and  N3  sites  independently.  Photoelectron  spec- 
troscopy'  "*-^  of  isolated  U“  ions  has  shown  that  the  extra  electron 
is  weaJcIy  bound  by  the  large  dipole  moment  of  U  (4.7  D).*^  In 
such  dipole-bound  states,  the  electron  is  localized  off  the  positive 
end  of  the  molecule  in  a  very  diffuse  orbital  consistent  with 
the  low  binding  energy.  In  a  different  experiment,^-^  electron 
transfer  from  Rydberg  state  Xe  to  U  and  U(Ar)„  clusters  was 
observed  for  electron  energies  of  0.01  “0.5  eV,  yielding  again 
the  dipole-bound  ions  U",  but  also  covalent- state  U"  as  well. 
Photoelectron  spectroscopy  of  cluster  anions  has  since  yielded 
an  estimate  of  the  binding  energy  of  the  covalently  bonded  extra 
electron  in  U“,  i.e.,  the  electron  affinity  of  U.^ 


♦  Corresponding  author.  E-mail:  thomas.miller@hanscom.af.mil. 
^  Air  Force  Research  Laboratory. 

*  ChemMotif,  Inc. 


The  initial  report  of  the  Rydberg  electron-transfer  study^ 
included  observation  of  (U-H)“  at  very  low  electron  energies, 
which  was  surprising  to  us  because  the  reaction 

e"  +  U  — (U-H)‘  +  H  (1) 

can  only  occur  (at  0  K,  at  least)  if  uracil  is  a  gas-phase 
superacid, that  is,  if  the  gas-phase  acidity  of  uracil  is  less 
than  Af/°acid  =  313.6  kcal  mol"*.  Another  way  of  saying  this 
is  to  view  reaction  I  as  a  proton-transfer  reaction  involving  the 
acids  U  and  H,  and  the  reaction  will  proceed  in  the  direction 
shown  only  if  U  is  a  stronger  gas-phase  acid  than  H.^*'  The 
possibility  of  dissociative  electron  attachment-  seemed  in  conflict 
with  the  observation  of  a  dipole-bound  parent  ion  because  the 
extra  electron  would  be  able  to  deprotonate  the  U  neutral  core 
via  reaction  1.  However,  we  could  not  rule  out  the  possibility 
that  the  most  acidic  site  of  U  might  not  be  at  the  positive  end 
of  the  molecule““there  are  four  hydrogen  atoms  to  consider^ 
or  that  a  large  conformational  barrier  exists  between  the  U" 
dipole-bound  state  and  deprotonation  of  U.  Calculations  indicate 
that  covalent  state  U“  is  a  puckered  ring  in  contrast  to  the  planar 
U.^  Although  the  initial  observation  of  reaction  1  was  with¬ 
drawn,^  we  were  intrigued  as  to  the  gas-phase  acidity  of 
nucleotide  bases  and  undertook  a  study  of  proton-transfer 
reactions  involving  U  and  of  electron  attachment  to  U  in  a  weak, 
thermally  equilibrated  plasma,  and  wc  carried  out  calculations 
of  bond  strengths  and  gas- phase  acidities  for  the  four  hydrogen 
sites  in  U.  Electron  attachment  to  nucleic  acid  bases  in  solution 
has  been  the  subject  of  a  number  of  studies  because  it  is  a  crucial 
step  in  the  evolution  of  radiation  damage  to  nucleic  acids.^^  In 
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TABLE  1:  Reactions  between  Reference  Acids  and  the 
Conjugate  Anion  of  Uracil,  Studied  with  the  SIFT 
Apparatus  at  298  K _ 

A//°acid''  (kcal  morO  rate  const^ 

ref  acid  AG^acw"  (kcal  mor‘)  reaction  channel  (10“^^  cm^  s~^) 


HF 

371,4  ±0.0 

adduct  ( l(X)%) 

<0.27 

365.5  ±  0,2 

coll  22.6 

HCl 

333.4  ±  0.0 

adduct  (75%) 

8.0 

328.1  ±0.2 

H'*'  transfer  (25%) 

coll  8.9 

HBr 

323.5  ±0.1 

H"^  transfer  (98%) 

7.8 

318.3  ±0.2 

adduct  (2%) 

coll  12.8 

°  Gas-phase  acidities  of  the  reference  acids  were  taken  from  ref  41. 
^  Estimated  accurate  to  within  25%,  The  calculated  collisional  rate 

constant  (see  text)  is  denoted  by 

“coll.” 

deoxyribonucleic  acid  (DNA),  measurements  indicate  that  most 
of  the  attachment  is  to  cytosine,  and  some  of  the  attachment  is 
to  thymine.  The  electron  moves  to  thymine  when  the  material 
is  annealed,  and  at  elevated  temperatures  the  thymine  anion 
protonates.^^  Presumably,  in  ribonucleic  acid,  uracil  plays  this 
same  role  in  radiation  damage.  There  are  a  number  of  studies 
on  the  electrical  nature  of  most  done  with  measure¬ 

ments  of  DNA  conductivity  in  solutions,  although  not  exclu- 
sively.^^’^® 

While  dipole-bound  anions  of  nucleic  acid  bases  have  been 
observed  in  the  gas  phase,  in  solution  covalent  or  valence 
anions  exist,  stabilized  by  the  binding  of  the  solvent  molecules. 
The  stabilizing  influence  of  ligands  has  been  demonstrated  in 
photoelectron  spectroscopy^’"^  and  Rydberg  state  atom  experi¬ 
ments.^ 

11.  Experimental  Section 

A.  SIFT.  A  selected-ion  flow  tube  (SIFT)  was  used  in  that 
portion  of  the  experiments  in  which  (U-H)“  ions  were  reacted 
with  reference  acids.  The  technique  and  this  apparatus  have  been 
amply  described  in  the  past.^*’^^  The  (U-H)“  ions  were  generated 
by  placing  uracil  powder  in  the  heated  Ar  reservoir  of  a  high- 
pressure  nozzle.  An  electron  beam  in  the  vacuum  just  outside 
the  exit  of  the  nozzle  source  created  a  plasma  of  mostly  Ar"*" 
and  electrons  in  the  gas  Jet,  from  which  (U-H)“  ions  were 
obtained  via  electron  attachment.  The  ions  were  mass  selected 
with  an  rf  quadrupole  mass  spectrometer  and  then  injected  into 
a  flow  tube  using  a  He  aspirator.  The  flow  tube  contained  a 
fast  flow  (125  m  s"^)  of  He  buffer  gas  at  53  Pa  pressure.  At 
either  of  two  inlets  (halfway  or  two-thirds  of  the  distance  down 
the  flow  tube),  a  known  flow  of  reference  acid  vapor  could  be 
added  to  the  flow  tube,  and  the  resulting  exponential  attenuation 
in  the  (U-H)“  intensity  was  measured  for  different  reference 
acid  concentrations.  The  flow  tube  was  terminated  with  an 
orifice  through  which  a  sample  of  the  ion  population  in  the  flow 
tube  was  passed  to  a  second  rf  quadrupole  mass  spectrometer 
at  high  vacuum.  Only  (U-H)“  ions  were  observed  in  the  mass 
spectra  when  no  reactant  gases  were  added  to  the  flow  tube. 
The  bulk  of  the  flow  tube  gases  were  pumped  away  with  a 
mechanical  pump.  Reference  acids  used  in  this  work  were  HF, 
HCl,  and  HBr.  The  SIFT  results  are  presented  in  Table  1. 

B.  FALP.  A  flowing-afterglow  Langmuir  probe  (FALP) 
apparatus  was  used  for  the  present  experiments  on  ion  reactions 
with  U.  This  apparatus^^’^"^  and  the  FALP  technique^^  have  been 
well  described  in  the  literature.  For  the  ion-molecule  reaction 
studies  reported  here,  the  FALP  was  used  as  a  conventional 
flowing  afterglow. In  this  case,  the  movable  Langmuir  probe 
played  a  limited  role:  to  make  sure  that  the  ion  density  was 
low  (<10^  cm“^)  enough  to  avoid  ion -“ion  recombination,  to 
make  sure  that  no  free  electrons  entered  the  reaction  zone,  and 
to  measure  the  ion  velocity.  The  electron  attachment  data  given 


here  utilize  the  full  capabilities  of  the  FALP,^^  in  which  the 
Langmuir  probe  is  used  to  measure  the  decay  of  the  electron 
density  as  attachment  events  occur. 

In  the  FALP  apparatus,  a  fast  flow  (ca.  1 1  000  std  cm^  min“^) 
of  helium  buffer  gas  was  established  in  a  flow  tube  of  3.7  cm 
radius.  All  data  below  were  obtained  with  a  flow  tube  pressure 
of  133  Pa  except  for  a  few  tests.  A  plasma  (c“,  He"^)  was  created 
at  the  upstream  end  of  the  flow  tube  with  a  microwave 
discharge.  Metastable-excited  helium  was  also  created,  but  it 
was  eliminated  through  Penning  ionization  by  addition  down¬ 
stream  of  Ar.  The  flow  velocity  ( 1 20  m  s"' )  of  the  plasma  was 
measured  using  the  Langmuir  probe  to  follow  the  propagation 
of  a  pulsed  disturbance  of  the  microwave  power  down  the  length 
of  the  flow  tube.  Uracil  vapor  was  added  at  a  point  about  57 
cm  downstream  of  the  microwave  discharge,  and  43  cm  from 
an  ion-sampling  aperture,  for  studies  of  both  ion  and  electron 
interactions  with  gas-phase  U.  Ions  exiling  the  flow  tube  through 
the  sampling  aperture  were  mass  analyzed  with  an  rf  quadrupole 
mass  filter  and  detected  with  a  continuous-channel  particle 
counter. 

Uracil  powder  (Aldrich,  98%)  was  loaded  into  three  glass 
tubes  (6  mm  outside  diameter)  running  along  the  length  of  the 
reaction  zone  against  the  wall  of  the  flow  tube.  The  ends  of  the 
glass  tubes  were  drawn  out  and  bent  along  a  radius  of  the  flow 
tube.  Thus,  the  uracil  inlet  consi.sted  of  three  needles  (1  mm 
inside  diameter)  protruding  toward  the  How  tube  axis  and 
terminating  about  I  cm  from  the  axis.  The  U  was  vaporized  by 
heating  the  entire  flow  tube  to  467  K,  at  which  point  the  vapor 
pressure  of  U  is  approximately  3.3  Pa.-^^  Helium  gas  was  passed 
through  the  glass  tubes  that  held  the  U  to  transport  the  vapor 
into  the  flow  tube.  With  no  helium  flowing  over  the  U,  the 
concentration  of  U  in  the  reaction  zone  was  about  10^®  cm”^. 
With  45  std  cm^  min~^  helium  flowing,  the  concentration  of  U 
rose  to  (1-2)  x  10^*  cm~^.  The  U  flow  could  be  reduced  to 
about  10^  cm“^  with  gentle  pumping  backward  through  the  glass 
tubes.  The  concentration  of  the  uracil  vapor  was  estimated  by 
adding  H2O  gas  upstream  and  reacting  the  resulting  H3O+  ions 
with  uracil  to  give  protonated  uracil,  U'H"^,  as  in  reaction  2. 

+  U  —  U-H+  +  H2O  +  44  kcal  tnoF'  (2) 

The  44  kcal  mol“*  exothermicity  of  the  reaction  is  the 
difference  in  the  proton  affinities  of  water  and  uracil.*^  If  the 
rate  constant  for  reaction  2  is  known,  then  measuring  the  HsO"^ 
decay  rate  can  be  used  to  determine  the  uracil  concentration  in 
the  flow  tube.  However,  since  the  rate  constant  for  reaction  2 
has  not  been  measured,  we  calculated  the  collisional  rate 
constant^^  based  on  the  gas-phase  dipole  moment*®  and  polar¬ 
izability^^  of  uracil.  There  exists  a  great  deal  of  data  showing 
that  exothermic  proton  transfer  reactions  usually  take  place  with 
a  rate  constant  that  is  at  or  near  the  collisional  value."**^  Because 
of  the  importance  of  H30'^  in  trace  gas  analysis,  it  is  particularly 
well  documented  the  that  the  rates  for  exothermic  proton  transfer 
between  HsO"^  and  organic  molecules  proceed  at  or  near  the 
collisonal  rates, so  this  method  should  provide  a  good  measure 
of  U  concentration  in  the  flow  tube. 

Ion~uracil  reactions  were  studied  by  introducing  a  source 
gas  into  the  afterglow  to  produce  negative  ions  that  were  the 
conjugate  bases  of  reference  acids.  Gases  used  for  this  purpose 
(and  in  square  brackets  the  negative  ions  produced  from  them) 
were  CF3CI  [Cl'],  CF3Br  [Br"],  CF3I  B”],  (CF3S02)0  [CF3S03~], 
CHF2CO2H  [CHF2CO2’],  cf3C(0x:h2C(0)CH3  [CF3C(0)chc- 
(0)CH3-],  CF3C(0)CH2C(0)CF3  [CF3C(0)CHC(0)CF3"],  and 
CF3C(0)SCH2CH3  [CF3C(0)S“].  In  cases  where  the  reference 
acid  A  itself  was  used  as  the  source  gas,  a  common  problem 
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TABLE  2:  Reactions  between  Uracil  and  Conjugate  Anions  of  Reference  Acids  (e.g.,  Cl"  from  HCl),  Studied  with  the  FALP 
Apparatus  at  467  K 


ref  ion  [acid] 

A//°acid"  (kcal  mol  ’) 
AG^acid**  (kcal  mor’) 

reaction  channel 

rate  const'' 

(10  ‘”cm^  s“') 

Cl  [HCi] 

333.4  ±  0.0 

transfer  (70%) 

23 

CHFzCOr  [CHF2CO2H] 

328.1  ±0.2 

adduct  (307r) 

coll  38.8 

331.0  ±2.2 

adduct  (55%) 

14 

CF,C(0)CHC(0)CH,-  [CFjC(0)CHC(OH)CH,] 

323.9  ±  2.0 

H+  transfer  (45%) 

coll  28.1 

328.4  ±  2.9 

transfer  (65%) 

3.1 

Br-  [HBr] 

321.9  ±2.0 

adduct  (35%) 

coll  25.0 

323.5  ±0.1 

adduct  (99%) 

23 

CFjQOS-  (CFjC(O)SHl 

318.3  ±0.2 

transfer  (1%) 

coll  29.5 

318.6  ±3.0 

adduct  (70%) 

1.0 

CF,C(0)CHC(0)CF,-  [CF,C(0)CHC(OH)CF,l 

312.6  ±2.0 

transfer  (30%) 

coll  26.0 

317.4  ±2.2 

310.2  ±2.0 

adduct  (100%) 

0.89 
coll  23.6 

1  [HI] 

314.4  ±0.0 

adduct  (95%) 

2.9 

CF5SO,-  [CF3S0,H] 

309.2  ±  0.2 

H+  transfer  (<  5%) 

coll  26.1 

305.4  ±  2.2 

adduct  (100%) 

0.81 

299.5  ±  2.0 

transfer  (<1%) 

coll  25.2 

“  Gas-phase  acidities  of  the  reference  acids  were  taken  from  ref  41.  "  Estimated  accurate  to  within  a  factor  of  2.  The  calculated  collisional  rate 
constant  (see  text)  is  denoted  by  “coll.” 


was  rapid  formation  of  a  cluster  anion  {A-H)“*A  in  addition  to 
the  desired  (A-H)".  In  these  cases,  the  acid  vapor  was  diluted 
with  Ar  gas  to  permit  more  accurate  adjustment  of  the  source 
gas  flow  rate  to  minimize  the  secondary  reaction  that  formed 
the  cluster  anion.  If  the  vapor  pressure  of  the  source  gas  was 
too  low  to  supply  a  concentration  in  the  flow  tube  adequate  to 
remove  all  free  electrons  prior  to  the  uracil  reaction  zone,  Ar 
gas  was  flowed  past  the  bulb  containing  the  source  gas  in  order 
to  carry  a  greater  amount  of  source  gas  into  the  flow  tube  by 
aspiration. 

High-resolution  mass  spectra  were  obtained  in  each  ion- 
uracil  reaction  case  to  identify  the  ionic  products  of  reaction. 
Low-resolution  mass  spectra  were  also  obtained  (to  minimize 
mass  discrimination  effects).  Peak  heights  in  the  low-resolution 
spectra  were  used  to  determine  the  reaction  rate  constant  using 
the  uracil  concentration  derived  from  the  H3O+  data  and  the 
reaction  time  measured  directly  as  described  earlier.  The  uracil 
concentration  and  reaction  time  are  merely  useful  intermediaries; 
the  reaction  rate  constants  reported  here  are  ultimately  normal¬ 
ized  to  the  amount  of  proton  transfer  observed  for  reaction  2 
for  the  same  U  concentration  and  reaction  time,  with  the 
assumption  that  the  rate  constant  for  reaction  2  is  collisional  in 
magnitude.  Since  the  rate  constant  for  reaction  2  has  not  yet 
been  measured,  we  do  not  claim  that  these  results  are  accurate 
better  than  a  factor  of  2.  However,  for  bracketing  the  gas-phase 
acidity  of  U,  we  are  mainly  concerned  with  two  questions:  does 
proton  transfer  from  uracil  to  the  reference  anion  take  place, 
and  is  the  reaction  fast  or  slow?  The  present  FALP  data  satisfy 
these  goals.  We  believe  that  the  relative  accuracy  of  the  reaction 
rate  constants  is  about  ±50%.  The  branching  fractions  for  ionic 
products  are  good  to  within  10  percentage  points.  The  results 
of  reaction  with  U  are  presented  in  Table  2. 

Secondary  clustering  reactions  were  noted  in  the  mass  spectra 
(and  accounted  for  in  the  product  analyses): 

(U-H)'"  +  U  +  He  —  (U-Hf-U  +  He  (3) 

(U-H)~  +  U  +  He  —  (U-H)”-U  +  He  (4) 

III.  Experimental  Gas  Phase  Acidity  Results 

Results  of  the  reaction  of  (U-H)"  with  HF,  HCl,  and  HBr  (at 
298  K)  are  given  in  Table  1.  Adduct  formation  occurs  in  all 
three  cases,  which  is  not  unexpected  since  adduct  formation  is 
enhanced  for  dipolar  molecules.  From  these  results,  it  is  clear 


that  HBr  is  a  stronger  acid  than  U  because  HBr  readily  donates 
a  proton  to  (U-H)“,  to  form  Br**  +  U.  It  is  also  clear  that  U  is 
a  weaker  acid  than  HF,  as  HF  will  not  donate  a  proton  to  (U- 
H)“.  The  HCl  case  is  in  the  middle  ground:  proton  transfer 
occurs,  but  at  a  relatively  slow  rate.  From  the  three  SIFT  data 
alone,  one  can  conclude  that  the  acidity  of  U  is  similar  to  that 
of  HCl,  for  which  A//°acid  is  333.4  ±  0.0  kcal  mol“‘,  and  AG'^acid 
is  328.1  ±  0.2  kcal  mol"*.^* 

Results  of  the  ion-uracil  reaction  experiments  (at  479  K) 
are  given  in  Table  2.  Proton  transfer  is  the  only  rearrangement 
channel  observed.  It  is  also  seen  from  Table  2  that  adduct 
formation  is  common,  even  dominant.  Data  were  obtained  with 
Cl“,  Br“,  and  I"  +  U  at  helium  buffer  pressures  of  57,  88, 
120,  and  133  Pa  to  see  if  we  were  in  the  high-pressure  limit  for 
stabilization  of  the  adduct  ion,  since  competition  from  the  adduct 
channel  complicates  an  otherwise  straightforward  interpretation 
of  the  results  in  Table  2.  We  found  that  the  branching  fraction 
for  the  adduct  channel  in  this  pressure  range  was  20—30%  for 
Cl“,  98—99%  for  Br“,  and  92—95%  for  1“.  Therefore,  while 
the  experiment  would  be  simpler  to  interpret  in  a  near-zero 
pressure  apparatus  (an  ion  cyclotron  resonance  cell  or  a  guided 
ion  beam  experiment),  operating  at  the  lowest  feasible  pressures 
in  an  FALP  was  not  very  helpful.  On  a  gross  scale,  it  is  clear 
in  Table  2  that  reference  acids  higher  on  the  list  readily 
deprotonate  uracil  while  those  lower  on  the  list  do  not. 

If  both  forward  and  reverse  reactions  are  observed  as  a  single 
channel,  then  the  gas-phase  acidity,  AC°acid»  can  be  calculated 
directly  from  the  ratio  of  the  forward  and  reverse  rate  constants. 
The  reaction  data  reported  in  Tables  1  and  2  are  complicated 
by  several  factors,  most  notably  that  the  forward  and  reverse 
reactions  were  studied  at  different  temperatures,  and  that 
clustering  reactions  compete  with  the  proton-transfer  reactions 
(see  product  ratios  in  Tables  1  and  2).  Therefore,  AG^acid  cannot 
be  calculated  from  the  ratio  of  forward  and  reverse  proton- 
transfer  rates,  even  when  they  are  observed  in  both  directions 
as  they  are  for  the  HCl/CI”  reference  acid^ase  pair. 

In  low-pressure  experiments,  the  kinetics  of  a  proton- transfer 
reaction  are  typically  governed  by  the  heal  of  deprotonation, 
A//°acid»  at  least  for  a  barrierless  reaction  such  as  this  simple 
proton  transfer  should  be.  However,  the  observation  of  the 
clustering  channel  suggests  that  the  observed  reactions  might 
not  meet  the  “low-pressure”  criterion  for  the  reaction  kinetics 
to  be  entirely  governed  by  A//°  rather  than  AG°.  That  is,  enough 
collisional  thermalization  by  the  buffer  gas  to  cause  clustering 
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AH°  .  of  reference  acid  (Real  mol’^) 


Figure  1.  Partial  reaction  rate  constants  for  proton  transfer  from  uracil 
(U)  to  conjugate  anions  (A-H)“  of  reference  acids  at  467  K  (solid 
circles)  and  for  proton  transfer  from  reference  acids  (A)  to  (U-H)“  at 
298  K  (solid  squares). 

reactions  also  means  that  the  pressure  is  approaching  that 
required  for  thermodynamic  equilibrium  conditions  to  hold.  Our 
past  experiences  with  reactions  for  which  association  competes 
with  a  chemical  channel  show  situations  in  which  clustering 
does  not  compete  with  proton  transfer^^  '^^  (therefore  A//°  would 
apply)  and  others  in  which  the  buffer  clearly  plays  a  role  in  the 
apparent  bi molecular  reaction  (and  therefore  AG°  would  ap¬ 
ply.)"^  It  was  not  possible  in  the  experiments  reported  here  to 
distinguish  whether  the  rates  for  the  proton-transfer  channel  are 
affected  by  collisions  with  the  buffer  gas. 

In  the  present  work  we  bracket  A//°acid(U)  from  the  proton- 
transfer  rates.  However,  the  magnitude  of  the  difference  between 
A//°acid  and  AG°acid  is  fairly  uniform  for  the  acids  used  in  the 
present  work,  and  the  slight  discrepancy  is  well  within  the 
accuracy  of  the  experiments.  Kinetics  of  the  proton-transfer 
reactions  observed  on  both  the  SIFT  and  FALP  apparatuses  are 
shown  in  Figure  1  as  a  function  of  A//°acid  of  the  reference 
species;  the  A//°acid  values  are  plotted  for  298  K.  Although  the 
acidity  for  each  species  is  slightly  larger  at  the  467  K 
temperature,  both  the  reference  and  U  acidities  shift  by  about 
the  same  amount:  a  0.6  kcal  mol"^  increase  in  A//°acid(UH) 
between  298  and  467  K  was  determined  from  Gaussian 
calculations  (discussed  below),  which  is  well  within  the  error 
limits  of  the  bracketing  experiments. 

The  results  in  Figure  1  can  be  said  to  bracket  A//°acid(kJ)  close 
to  that  of  HCl,  or  333  kcal  mol'^  Using  the  entropy  correction 
calculated  in  section  IV,  this  result  implies  AG°acid(U)  =  326 
±  5  kcal  mol“^ 

The  relation  between  the  gas-phase  A//°acid  and  the  solution- 
phase  acidity,  or  pJ^a,  may  be  obtained  from  a  thermochemical 
cycle  and  is  given  by  eq  5. 

=  23RT(pK,)  +  TAS\  +  AH\^(U)  - 

A//V[(U-H)“]  -  AH\^iH^)  (5) 

The  term  23RT(pKa)  gives  the  AG°  for  the  acidity  in  aqueous 
solution  (pKa  =  9.5;"^^  AG°  =  +12.95  kcal  mol"^.  The  7AS°a 
term  is  the  entropy  for  aqueous  acidity;  an  entropy  for  the 
solution-phase  acidity  of  20  cal  mol“^  is  assumed, yielding 
a  net  AH°^  in  solution  of  +18.9  kcal  mol“^  The  AH\y^  terms 
are  the  enthalpies  of  hydration  of  the  various  species,  i.e.,  the 
energetics  for  transfer  from  gas  phase  to  solution.  A//°hydr  (U) 


is  taken  as  the  sum  of  the  aqueous  solvation  enthalpy  of  solid 
uracil,  +7.0  ±  0.3  kcal  mo^^'*^  and  the  negative  of  the 
sublimation  enthalpy  of  solid  uracil,  -31.3  ±  1.2  kcal  moU^"^”^ 
for  a  net  hydration  enthalpy  of  —24.3  kcal  mol“^  The  hydration 
enthalpy  for  the  proton  is  taken  to  be  —262  kcal  mol"*. 
Substitution  of  the  known  values  into  eq  4  yields  a  solvation 
enthalpy,  A//°hydr[(U-H)"],  of  —76.4  kcal  mol"*. 

The  hydration  enthalpy  of  a  spherical  charge  is  given  by  the 
Bom  theory  of  solvation,  with  the  result  in  eq  6  of 

A//“hydr  =  -  (1  -  = 

—  \65.7/r^  kcal  mol”’  (for  in  A)  (6) 

in  which  M  is  the  number  of  ions  solvated,  qi  is  the  charge  on 
the  ions,  L  is  a  small,  temperature-dependent  correction  to  the 
dielectric  constant  (4.63  x  10"^  for  water).  T  is  the  temperature, 
D  is  the  dielectric  constant  (78,53  for  water),  and  n  is  the 
solvation  radius  for  the  ion.^*  Substituting  AA/°hydr[(U-H)"]  = 
—76.4  kcal  mol"*  into  eq  5  gives  a  2. 1  A  solvation  radius  for 
(U-H)"  in  aqueous  solution.  This  is  a  surprisingly  small  number 
for  the  effective  solvation  radius,  considering  both  the  size  of 
(U-H)"  (see  section  IV)  and  the  fact  that,  for  anions,  the 
solvation  radius  is  typically  about  0.85  A  larger  than  the  crystal 
radius  of  the  ion.^*  This  suggests  that  the  hydration  of  (U-H)“ 
is  helped  considerably  by  other  elTecls,  most  probably  very 
effective  hydrogen  bonding  to  the  oxygens,  to  the  anionic  Nl, 
and  from  the  N3  hydrogen.  The  hydration  enthalpy  for  the  parent 
U  is  somewhat  large  as  well,  by  comparison  to  other  similarly 
sized  organic  molecules,  and  is  again  indicative  of  strong 
interactions  between  the  U  and  the  H2O  of  solvation. 

No  U"  was  observed  in  either  the  SIFT  or  FALP  work,  that 
is,  within  the  uncertainty  of  measurement  of  deviation  from  the 
expected  *^C  component  of  U"  at  112  amu.  Even  if  covalent- 
state  U"  were  produced,  we  would  not  expect  a  polyatomic 
negative  ion  of  such  low  electron  affinity  (3.5  ±  2.8  kcal  mol"*)^ 
to  survive  in  the  buffer  gas  at  the  pressures  and  temperatures 
of  the  present  flow  tube  experiments.  The  dipole-bound  U", 
with  its  lower  binding  energy  (2.0  ±  0.2  kcal  mol"*),^  would 
stand  still  less  chance  of  surviving. 

The  FALP  apparatus  was  also  used  to  investigate  reaction  1. 
The  Langmuir  probe  was  used  to  measure  the  decay  in  the 
electron  density  along  the  axis  of  the  flow  tube  when  U  vapor 
was  present.  Very  few  negative  ions  were  produced  by  attach¬ 
ment.  While  (U-H)"  ions  were  observed,  there  was  little  decay 
in  the  electron  density  along  the  flow  tube  aside  from  that  due 
to  ambipolar  diffusion.  An  upper  limit  of  2  x  10"*®  cm^  s"* 
may  be  placed  on  the  electron  attachment  rate  constant  at  467 
K.  The  very  weak  electron  attachment  signal,  even  if  partly  due 
to  impurities,  is  consistent  with  a  large  endothermicity  (17  kcal 
mol"*  according  to  this  work)  for  reaction  1.  Experiments  were 
also  carried  out  with  suprathermal  electrons.  For  this  purpose, 
pure  At  was  used  as  the  buffer  gas,  in  which  case  the  electron 
temperature  is  about  2000  K  because  Ar  does  not  thermalize 
low-energy  electrons."*^  These  electrons  have  only  a  small 
fraction  of  the  energy  needed  to  overcome  the  endothermicity 
of  reaction  1,  so  it  was  not  surprising  that  electron  attachment 
remained  an  immeasurably  weak  process.  This  observation  may 
be  compared  with  results  of  electron  beam  experiments  where 
major  resonances  were  observed  only  at  energies  above  1  eV."*^ 
A  smaller  feature  for  (U-H)"  production  was  observed  around 
0.65  eV,  which  can  be  explained  by  internal  energy  in  the  target 
molecules.  However,  it  is  difficult  to  rationalize  a  weak  feature 
observed  at  about  60  meV  as  being  due  to  (U-H)".  A  large 
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Figure  2.  Structures  of  uracil  molecule  (lop),  neutral  (U-H)  (bottom  left,  H  atom  missing  from  the  most  acidic  site),  and  (U-H)“  anion  (bottom 
right,  missing  from  the  most  acidic  site).  The  arrow  shows  the  direction  of  the  dipole  moment  (4.7  D)  in  U.  The  C-C,  C-N,  and  C-0  bond 
lengths  were  determined  from  MP2(Full)/6-31 1++G(3df,2p)  calculations;  those  given  in  italics  are  the  MP2(Fuil)/6-31G(d)  geometrv  used  in  the 
G3  and  G3(MP2)  calculations. 


zero-energy  resonance  was  originally  reported  but  later  found 
to  be  due  to  the  calibration  gas  used  in  the  experiment. 

IV.  Gaussian  Calculations 

Electronic  structure  calculations  were  carried  out  with  the 
primary  aim  of  determining  which  of  the  four  H-atom  sites  in 
uracil  is  most  acidic  in  the  gas  phase.  The  calculations  include 
the  homolytic  bond  strengths  and  the  electron  affinities  of  the 
four  possible  U  radicals,  information  which  aids  in  understand¬ 
ing  the  acidic  strength.  G3  and  G3(MP2)  calculations  were 
performed  on  the  structures  and  energies  of  U,  (U-H),  and  (U- 
H)”  using  the  GAUSSIAN-98W  package  of  programs. Fol¬ 
lowing  the  G3  and  G3(MP2)  formalisms,^-’^^  molecular  struc¬ 
tures  were  first  optimized  at  the  HF/6-3IG(d)  level  of  theory, 
and  the  stability  of  the  wave  functions  were  checked  (i.e.,  the 
electron  configuration  was  altered  in  various  ways  to  verify  that 
the  lowest  energy  state  had  been  located).  Then  vibrational 
energies  were  calculated,  scaled  by  an  empirical  factor  of 
0.8929,^^  and  from  these  and  rotational  energies  the  zero-point 
energies  were  determined  along  with  thermal  corrections  to 
energies  (with  the  standard  state  taken  as  298.15  K  and  1  atm). 
Next,  structures  were  optimized  at  the  MP2(FuIl)/6-3  lG(d)  level. 
The  structure  of  U  was  also  optimized  with  a  large  basis  set  [at 
the  MP2(Full)/6-31 1++G(3df,2p)  level]  for  comparison  with 
the  results  obtained  with  the  smaller  basis  set  used  in  G3  theory, 
and  with  the  moderate  6-31  H"4'G(d,p)  basis  set  used  by 
Dolgounitcheva  et  al.*^  The  use  of  many  polarization  functions 
in  the  large  basis  set  tended  to  shorten  bond  lengths  by  as  much 
as  O.Ol  A  compared  to  those  obtained  with  the  6-31G(d)  basis 
set.  The  bond  lengths  are  given  in  Figure  2,  which  shows  the 
structures  of  U  and  (U-H)“  relevant  to  the  acidity  work. 
Cartesian  coordinates  for  U  and  all  (U-H)  and  (U-H)"  are  given 
in  Supporting  Information.  Uracil  is  known  to  be  planar  and  to 
have  the  diketo  tautomer  structure  in  the  gas  phase.^*^^  The 
geometry  optimizations  were  initially  restricted  to  planar  U, 


(U-H),  and  (U-H)".  Small  out-of-plane  distortions  were  then 
introduced  (typically  10°),  and  MP2  optimization  of  the  distorted 
geometries  led  back  to  the  planar  forms.  The  one  exception  is 
neutral  (U-H)  with  an  H  atom  missing  from  the  N3  site.  The 
optimizations  and  frequency  analyses  consistently  led  to  a 
slightly  nonplanar  geometry,  with  appendage  atoms  pro¬ 
truding  from  1°  to  5°  above  and  below  the  average  plane  of  the 
ring.  G3  calculations  showed  the  nonplanar  form  to  be  0.13 
kcal  mol"*  lower  in  energy  than  the  planar  one  (^A').  DFT 
optimization,  with  a  larger  basis  set  than  used  in  G3,  led  to  the 
planar  form.  Thus,  a  small  academic  uncertainty  remains  at  the 
N3  site.  There  was  little  change  in  the  structures  calculated  for 
the  uracil  molecule  with  and  without  electron  correlation,  and 
increasing  account  of  correlation,  presumably  because  of  the 
fundamental  rigidity  of  a  six-membered  ring. 

Gas  phase  acidity  values  were  calculated  as  follows:  (a)  at 
0  K,  A//°acid  and  AG°acid  are  the  same,  simply  the  sum  of  the 
total  energies  of  the  separated  fragments  (U-H)“  +  H"^,  minus 
the  total  energy  of  U  itself;  (b)  A//°acid  at  298  K  is  the  difference 
in  the  enthalpy  of  U  and  that  of  the  separated  fragments;  and 
(c)  AG°acid  at  298  K  is  the  difference  in  the  free  energy  of  U 
and  that  of  the  separated  fragments.  The  total  energy  of  the 
proton  is  zero  at  0  K,  and  its  enthalpy  is  5kT12  at  298  K.  The 
homolytic  bond  su-engths  were  calculated  similarly,  for  separa¬ 
tion  of  U  into  neutral  fragments  (U-H)  +  H.  The  total  energy 
of  H  atom  is  0.5  hartree  at  0  K,  while  the  enthalpy  at  298  K 
includes  an  additional  5kT/2. 

Calculations  based  on  G3  and  G3(MP2)  theory  contain  two 
empirical  corrections:  one  is  the  scaling  of  Hartree -Fock  (HF) 
zero-point  corrections,  and  a  second  is  a  “higher  level  correc¬ 
tion”  which  accounts  for  remaining  deficiencies  in  the  methods 
through  a  comparison  of  calculated  and  measured  ionization 
potentials  and  electron  affinities.  The  higher  level  correction 
cancels  out  (exactly)  when  acidities  are  calculated.^^*-''^  Similarly, 
any  systematic  error  in  the  zero  point  energy  (ZPE)  results  tends 
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TABLE  3:  G3  and  G3(MP2)  Energies  and  Zero-Point  Energies  (ZPE),  in  hartrees,  for  Acidic  Sites  in  Uracil 
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G3(MP2) 


molecule^ 

ZPE* 

total  energy  (0  K) 

enthalpy  (298  K) 

free  energy  (298  K) 

total  energy  (0  K)^ 

uracil  (C4H4N2O2,  'A') 

0.084  487 

-414.571  427 

-414.564  168 

-414.601  919 

-414.257  772 

uracil-H,  N1  site  (C4H3N2O2,  ^A") 

0.069  697 

-413.911  554 

-413.904  266 

-413.942  875 

-413.596  729 

uracil-H+,  N1  site  (C4H3N202",  'A') 

0.071  509 

-414.040  397 

-414.033  500 

-414.070  579 

-413.726  714 

uracil-H,  N3  site^  (C4H3N2O2,  W') 

0.071  376 

-413.869  643 

-413.862  612 

-413,900  564 

-413.555  982 

uracil-H+,  N3  site  (C4H3N2O2-,  ‘A') 

0.070  801 

-414.021  363 

-414.014  146 

-414.051  994 

-413.707  812 

uracil-H,  C5  site  (C4H3N2O2,  ^A') 

0.071  575 

-413.881  486 

-413.874  162 

-413.912  695 

-413.566  783 

uracil-H+,  C5  site  (C4H3N2O2-,  'AO 

0.070  841 

-413,968  540 

-413,961  143 

-413,999  480 

-413.655  764 

uracil-H,  C6  site  (C4H3N2O2,  W) 

0.071  344 

-413.891  238 

-413.883  920 

-413,922  408 

-413.576  651 

uracil-H''’,  C6  site  (C4H3N202',  'AO 

0.071  428 

-413.990  114 

-413.983  031 

-414.020  433 

-413.676  895 

The  label  “uracil-H”  means  one  H  atom  removed  from  uracil;  “uracil-H+”  means  one  proton  removed  from  uracil.  *  ZPEs  and  thermal  corrections 
based  on  HF/6-31G(d)  frequencies  scaled  by  0.8929.  ^  G3(MP2)  enthalpies  and  free  energies  may  be  obtained  from  the  0  K  energy  using  differences 
in  the  G3  entries;  see  text.  ^  See  text:  the  nonplanar  form  differs  from  the  planar  ^A'  form  by  only  0.13  kcal  mol"^  and  is  therefore  uncertain. 


TABLE  4:  Calculated  Homolytic  Bond  Enthalpy  D°298[G4H3N202-H],  Electron  Affinity  (EA)  of  the  Radical  C4H3N2O2,  and 
Gas-Phase  Acidity  of  Uracil,  in  kcal  mol“^ _ 


quantity 

G3(MP2)  (0  K) 

G3{MP2)  (298  K) 

G3  (0  K) 

G3  (298  K) 

G3  (467  K)“ 

at  the  N1  vertex 

homolytic  bond  enthalpy 
EA(C4H3N202) 

101.1 

81.6 

102.6 

100.3 

80.9 

101.8 

102.5 

A//°acid(uracil) 

333.2 

334.5 

333.2 

334.5 

335.1 

AC°acid(uracil) 

333.2 

327.2 

333.2 

327.1 

327.6 

at  the  N3  vertex 

homolytic  bond  enthalpy 
EA(C4H3N202) 

126.6 

95.3 

128.0 

126.6 

95.2 

128.0 

A//°acid(uracil) 

345.1 

346.6 

345.1 

346.6 

AG°acid(uracil) 

345.1 

338.7 

345.1 

338.7 

at  the  C5  vertex 

homolytic  bond  enthalpy 
EA(C4H3N202) 

119.8 

55,8 

121.4 

119.2 

54.6 

120.7 

A//°acid(uracil) 

377.8 

379.3 

378.3 

379.9 

AG°acid(uracil) 

377,8 

371.2 

378.3 

371.8 

at  the  C6  vertex 

homolytic  bond  enthalpy 
EA(C4H3N202) 

113.7 

62.9 

115.2 

113.1 

62.0 

114.6 

A//°acid(uracil) 

364,5 

365.9 

364.8 

366.1 

AG°acid(uracil) 

364.5 

358.3 

364.8 

358.6 

°  Values  at  467  K  are  needed  for  interpretation  of  the  FALP  data;  see  text. 


to  cancel  when  differences  in  total  energies  are  taken  in  order 
to  obtain  acidities,  EAs,  and  bond  strengths.  Smith  and  Radom^^ 
and  Notario  et  al.^^  have  shown  that  gas-phase  acidities 
calculated  at  the  G2  level^^’^®  are  accurate  typically  within  ±2 
kcal  mol”^;  G3  theory  should  yield  an  improvement  over  this, 
approaching  ±1  kcal  mol"^  on  average,  G3(MP2)  error  should 
lie  between  those  of  G3  and  G2. 

The  G3  and  G3(MP2)  ZPEs  and  energies  are  given  in  Table 
3.  For  G3(MP2)  only  the  0  K  total  energy  is  listed,  but  the 
G3(MP2)  298  K  enthalpy  and  free  energy  may  be  obtained  from 
differences  in  the  G3  values  since  the  same  thermal  corrections 
are  used  in  G3  and  G3(MP2)  theory.  To  give  one  example,  for 
U,  the  G3(MP2)  total  energy  at  0  K  is  —414.257  77  hartrees. 
The  G3  numbers  imply  that  the  enthalpy  at  298  K  lies  0.007  26 
hartree  higher.  The  G3(MP2)  enthalpy  is  equivalently  higher, 
and  is  thus  —414.250  51  hartrees.  Gas-phase  acidities,  EA(U), 
and  homolytic  bond  enthalpies,  derived  from  the  energies  in 
Table  3,  are  given  in  Table  4.  Values  at  467  K  were  included 
for  interpretation  of  the  FALP  data.  The  calculated  acidity  (at 
the  N1  site)  is  in  agreement  with  the  experimental  value 
determined  in  this  work. 

The  results  at  0  K  necessarily  obey  the  thermochemical  cycle 
given  in  eq  7 

K)  =  D°o[(U-H)-H]  -  EA[(U-H)]  +  IE(H) 

(7) 


where  IE(H)  is  the  ionization  energy  of  H  atom.  Applying  this 
cycle  to  the  bond  enthalpies  and  bond  dissociation  enthalpies, 
with  all  quantities  given  in  Table  4,  aids  in  understanding  the 
moderate  acidity  of  uracil:  the  homolytic  bond  strength  is  not 
especially  weak,  but  the  EA  of  the  fragment  is  high. 

The  calculations  indicate  that  the  (U-H)”  ion  observed  in  this 
experiment  results  from  loss  of  from  the  N1  vertex  in  Figure 
2.  The  acidities  resulting  from  loss  of  alternative  protons  are 
too  large  to  play  a  role  in  the  thermal  energies  of  the  present 
experiments.  The  large  EA  calculated  for  (U-H),  81.9  kcal  mol“^ 
at  the  N1  site,  agrees  with  an  experimental  observation^  EA- 
[(U-H)]  >  57.6  kcal  mor^ 

G3  theory  is  intended  to  yield  a  good  approximation  to  a 
full  configuration  interaction  calculation.  G3(MP2)  calculations 
were  carried  out  for  comparison  with  the  G3  results  because 
the  G3(MP2)  ones  are  far  more  economical.  In  the  present  woik, 
G3(MP2)  calculations  required  1/7  to  1/5  the  computational  time 
of  the  G3  ones.  The  results  in  Table  4  show  good  agreement 
between  the  two  methods,  within  0.4  kcal  mol”^  on  average. 
The  calculations  described  here  were  initially  carried  out  at  the 
G2(MP2)  and  G2  levels^^^ss  theory.  The  G2(MP2)  and  G2 
results  (given  in  the  Supporting  Information)  agree  with  the 
corresponding  G3(MP2)  and  G3  results  within  0.5  kcal  mol"* 
on  average,  with  the  worst  discrepancy  being  EA(U-H)  at  the 
N3  site,  for  which  the  G2(MP2)  result  was  1.4  kcal  mol"^  lower 
than  the  G3(MP2)  value. 


Acidity  of  Uracil 
V.  Conclusions 

We  have  measured  rate  constants  for  proton-transfer  reactions 
between  the  conjugate  anions  of  reference  acids  and  uracil  (U), 
and  between  (U-H)“  and  reference  acids,  in  the  gas  phase.  The 
results  reported  here  show  that  the  gas-phase  acidity  of  uracil 
is  A//®acid(U)  =  333  ±  5  kcal  mol“*  [and,  by  implication, 
AC°acid(U)  =  326  ±  5  kcal  mor^]  at  298  K.  This  result  is  in 
complete  agreement  with  that  of  Kurinovich  and  Lee,^  who  used 
a  Fourier  transform  mass  spectrometer  (FTMS)  to  determine 
when  proton-transfer  reactions  were  taking  place  in  ionized  gas 
mixtures.  While  they  did  not  measure  reaction  rate  constants, 
their  bracketing  was  more  straightforward  than  in  our  experi¬ 
ments  because  the  gas  pressure  in  the  FTMS  is  low  enough 
that  there  is  no  competition  from  collisional  stabilization  of  the 
ion-molecule  complex.  They  measured  the  acidities  of  the  N1 
site.  Af/°acid(U)  =  333  ±  4  kcal  mol”',  and  the  N3  site,  347  ± 
4  kcal  mor’.  They  used  density  functional  methods  to  calculate 
acidities  of  376  and  362  kcal  mol"*  at  the  C5  and  C6  sites  in 
the  ring.  The  G3  calculations  in  the  present  work  should  be 
about  10  times  more  accurate  on  average  than  the  density 
functional  ones.  They  confirm  that  the  most  acidic  site  is  indeed 
at  the  N1  site  in  the  ring,  and  yield  A//°acid(U)  —  334.5  kcal 
mol”'  and  AG®acid(U)  =  327.1  kcal  mol”',  in  good  agreement 
with  the  present  experimental  result  and  that  of  Kurinovich  and 
Lee.^  The  G3  acidities  for  the  other  sites  are  A//°acid  =  346.6 
kcal  mo!”'  at  the  N3  site  (in  excellent  agreement  with 
Kurinovich  and  Lee^),  379.9  kcal  mol”'  at  the  C5  site,  and  366. 1 
kcal  mol”'  at  the  C6  site.  The  N1  hydrogen  is  of  course  lost  in 
RNA  formation,  since  the  N1  site  becomes  covalently  bonded 
to  one  of  the  carbons  of  the  ribose  sugar;  the  N3  site  is  hydrogen 
bonded  to  adenine.  The  G3  calculations  also  yield  homolytic 
bond  enthalpies  for  the  several  hydrogen  sites:  101 .8  kcal  mol"' 
(N1  site),  128,0  kcal  mol”'  (N3  site),  120.7  kcal  mol”'  (C5 
site),  and  1 14.6  kcal  mol”'  (C6  site). 

Our  calculated  G3  bond  energies  and  electron  affinities  for 
each  of  the  H-atom  sites  compare  quite  well  with  those  of  Hanel 
et  al.,^^  who  used  a  slightly  lower  level  of  theory  [G2(MP2)]. 

Electron  attachment  to  U  at  467  K  was  immeasurably  weak, 
consistent  with  the  endothermicity  implied  by  the  value  of  the 
gas-phase  acidity  determined  here. 
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